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ABSTRACT 



In anticipation of the AQN 201 nuclear reactor 0 s operation at a 
power of 1,000 watts for short periods of time and at 20 watts continuous- 
ly, a prediction of the behaviour during a rapid power rise is obtained 
from digital computer solutions of the reactor kinetic equations. The 
experimental and analytical work which provided the coefficients for the 
equations is outlined. Solutions of the equations indicate the peak power 
which may be achieved, the time during which 1,000 watts may be obtained, 
and the temperature rise that results. 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by his thesis advisor Associate Professor H. E. 
Handler, and by Associate Professor C. P. Howard and Associate Professor 
R. M. Thatcher, all of the U. S. Naval Postgraduate School, and for the 
assistance of the reactor facility personnel in obtaining experimental 
data. 
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1. Introduct ion. 



The AGN 201 nuclear reactor Is a low power reactor designed for 

educational uses. Usually, it is licensed to operate at 100 milliwatts 

and serves to demonstrate many of the basic characteristics of reactors® 

Also, at this power, the peak neutron flux of about 4,5 x 10 neutrons/ 

2 

cm -sec is sufficient for the creation of many radioactive isotopes in 
practical quantities® 

For other isotopes, however, it is desirable or necessary to have a 
much higher neutron flux and it has been proposed that the AGN 201 be 
operated for short periods at a power of 1,000 watts® Several of these 
reactors have been operated at somewhat higher powers, but none have 
operated as high as 1,000 watts. The limiting factor on elevated power 
operation is the lack of a cooling system. With no means of increasing 
the rate of heat transfer from the core, the overall core temperature in- 
creases. Since the reactor has a negative temperature coefficient of 
reactivity, a rise in temperature reduces the available excess reactivity, 
and the rate at which power may be further increased is reduced. 

The maximum steady- state power is that power at which the rate of 
heat generation in the core is equal to the rate of heat transfer from 
the core and the corresponding overall temperature of the core has re- 
duced the initial excess reactivity to zero. For the AGN 201, maximum 
steady-state power has been estimated to be approximately 20 watts. 

High powers, such as 1,000 watts, would be achieved during the transient 
period following a relatively large increase in reactivity if the re- 
actor core is initially at room temperature. Under such conditions, pow^r 
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can be increased rapidly by several orders of magnitude before the 
temperature has risen sufficiently to make the reactor sub-critical* 

The transient behaviour of a reactor can be predicted using reactor 
kinetic equations based on various approximations 0 In general, the equa- 
tions cannot be solved analytically due to the coupling between them, 
and analog computer solutions are frequently obtained for cases where 
variations of the dependent variables are not large* When the dependent 
variables change by orders of magnitude, numerical methods may be ap- 
plied to solve the equations using a digital computer,. 

It is the purpose of this investigation to establish a set of re- 
actor kinetic equations for the AGN 201 and obtain a solution to the 
equations using a digital computer for the case of a rapid rise in 
power. Features of significance to be obtained from the solution in- 
clude an estimate of the peak power that may be achieved (or in the run- 
away accident sense, a peak power which will not be exceeded) , the length 
of time that a power of 1,000 watts may be maintained, and the tempera- 
tures which would result. The latter is significant since there is a 
thermal fuse in the core which allows the bottom half of the core to drop 
away in the event of overheating. The fuse is always at a higher tempera- 
ture than the rest of the core. 

The coefficients of the kinetic equations are obtained from experiment- 
al data or from calculations. The available excess reactivity of the re- 
actor, the rate at which reactivity can be added, and an estimate of the 
temperature coefficient of reactivity can be obtained experimentally. 

The remaining nuclear coefficients and heat transfer coefficients must 
be calculated from the tabulated data for materials and compared where 
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possible with known physical characteristics of the reactor to estab- 
lish their reliability a For example* the nuclear parameters should 
predict the actual critical mass* and the calculated kinetic equa- 
tion coefficients* when applied to a computer solution for low power 
transient behaviour* should provide curves in reasonable agreement with 
curves obtained experimental ly. 

The dependent variables in the kinetic equations are temperature* 

neutron density* and the concentrations of the six delayed neutron 

precursors. The effect of nuclear poisons is omitted since the peak 

11 2 

flux at 1*000 watts is in the order of 10 neutrons/cm -sec at which 
poison effects are not significant. From neutron density* neutron flux 
and power can be determined. In the following sections* a description 
of the reactor* the procedure for obtaining parameters* and the method 
of obtaining computer solutions for neutron density and temperature 
described. 
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2. Description of the Reactor 

The reactor assembly is shown in Figure 1 and further detail of 
the core is shown in Figure 2. The core material is a homogeneous 
mixture of polyethylene and uranium dioxide; the uranium is 20 per- 
cent enriched in Uranium 235* The core is surrounded by a high- 
density graphite reflector, a lead shield, and a water shield* For 
powers higher than 100 milliwatts additional shielding is added in 
the form of concrete blocks* The structural material in the core 
region is 6061-T6 aluminum* 

Control is affected by the insertion and withdrawal of rods con- 
taining essentially the same material as the core. Three rods are 
the same size; two of these serve as safety rods, and the third is 
a coarse control rod* The fourth rod, which is smaller, provides 
fine control* 

Table 1 summarizes the characteristics of the AGN 201 series 
of reactors in general* 
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Table 1 AGN Characteristics 
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Figure 1. AGN 201 Reactor Assembly 
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Figure 2 0 AGN 201 Core Assembly 



3 a Preliminary Considerations* 

From the characteristics noted in the previous section* it is evident 
that the reactor has little excess reactivity at room temperature*, Dur- 
ing laboratory periods* it has been found that five or six of the small 
cadmium covers used around foils in thermal flux determinations are suf- 
ficient to shut the reactor down*, This limited excess reactivity and the 
strong negative temperature coefficient would prevent continuous operation 
with the glory hole empty at core temperatures above 27°C 0 Since the re- 
actor has no cooling system* the average temperature in the core would in- 
crease seven degrees quite rapidly when the reactor is operating with a 
high neutron flux and the resulting high heat generation rate in the 
core. 

One would expect then that it is necessary to approach high powers 
as rapidly as possible* that temperature would increase significantly 
with a corresponding reduction in reactivity when the reactor is at suf- 
ficiently elevated powers* and that in a relatively short time a negative 
reactivity would result causing a decrease in power 0 Subsequent long 
period oscillations of power about an equilibrium position probably would 
be the end result 0 

To establish a set of equations which predict this behaviour and 
can be solved* it is necessary to use a number of approximations* Fore* 
most of these is the Diffusion Equation approximation of the exact 
Boltzman Equation; in this approximation isotropy of the neutron veloci- 
ties is assumed o Also* exact energy distributions are replaced by the 
appropriate averages over energy intervals* Dividing the neutron popu- 
lation into groups according to the energy provides excellent results ii 
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a large number of groups are established. The complexity of solving 
this problem*, however, is not justified except in the most refined deter- 
minations of neutron behaviour in materials. 

One-group theory, modified one-group theory, Fermi Age theory and 
two-group theory are frequently used as approximations. Of these, Fermi 
Age theory is a good approximation for reactors moderated by relatively 
heavy materials, while two-group theory provides reasonable results for 
reactors in which hydrogen is the principal moderating material. (1) 

When the moderator contains hydrogen and a number of heavier materials, 
the Selengut - Goertzel equation yields better results than a two-group 
approach since the Selengut - Goertzel equation treats hydrogen slowing 
down and slowing down due to other materials separately. However, to 
obtain reasonable parameters for the time dependent problem with more 
than one region, two-group theory is more readily applied. 

A second area of approximations is the geometric model of the re- 
actor to be assumed. Since the theories require the solution of several 
simultaneous partial differential equations containing the Laplacian and 
two or three independent variables, a solution for the actual geometry 
would be extremely difficult to obtain. The problem is sufficiently 
complicated if one assumes a core region surrounded by a reflector re- 
gion with cylindrical geometry, and ignores all details of structural 
materials and construction. This model is satisfactory from the nuclear 
standpoint, but further consideration is required in connection with heat 
transfer, since the structural aluminum is a good conductor and the 
boundary and interface conditions have a different physical significance. 
These features are discussed in Section 6. 
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For a single region containing one material , the t ime- dependent 
neutron diffusion equation for a single energy group simplifies to: 

DV^Rt) -E<kf, t) + SP.t) = 

and for this same region the time-dependent heat diffusion equation is 
given as: 

KV 2 T(?,t) + =/°° c ^^ ] 



In both equations, the coefficients may be functions of the independent 
variables in the two equations; for example, the conductivity and dif- 
fusion coefficient are both temperature dependent • The source term in 
the heat diffusion equation is proportional to the thermal neutron flux 
for the same point in the region, and, in turn, the source term and the 
coefficients of the neutron diffusion equation are affected by the tempera- 
ture o 



Considering a small volume element in the region using a two-group 
approach, we have "the following coupled partial differential equations: 

D f v\ - E A + V EA(i -B) + Ct = I? (I = 1 to 6) 

_ b n c 



D s v 2 i - z a k * P IA = 

v eaa- x,.c v =-t£ 



KV 2 T +■ E<i> s =A c If 

where the subscripts "I" and n s" refer to the fast and slow neutron 
groups, respectively, and where E is a factor to convert neutron flux 
in the element to the corresponding heat generation rate 0 These equa- 
tions could be solved by numerical methods using finite difference 
techniques and a mesh system 0 The complexity of such a solution and 
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the time required to solve this problem are prohibitive even with the use 
of a high speed digital computer* Further discussion of this problem is 
included in Section 8, 

A more practical approach, and the approach used most frequently, is 
to obtain from the above equations a set of ordinary differential equa- 
tions amenable to computer solutions* Considering only the thermal group 
of neutrons, the diffusion equation may be written* 0 



non-leakage probability of the delayed neutrons* With the assumption 
that the flux can be calculated from the wave equation 



This equation coupled with the six ordinary differential equations for 
the delayed neutron precursors and a differential equation for tempera 
ture can be solved readily by numerical methods* 

Irrespective of the approach, it is necessary to obtain values for 
coefficients and parameters which are required for the equations* Some 
of these parameters can be calculated quite accurately* and some data 
are tabulated* In other cases, the values must be estimated* The test 
of any set of parameters selected is that they predict features of the 
reactor which are known* For example, the parameters should, when used 







the equation above may be reduced to 




(See Appendix 1) 
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in the appropriate equations, yield approximately the known critical mass 
and the steady-state flux distribution, and solutions of the transient 
equations at low power should agree with corresponding curves from the 
reactor. For this reason a number of experiments were carried out on 
the reactor to obtain data both to supplement calculated data and to 
provide verification of calculated data. These experiments are described 
in the following section* 



t 
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4* Experimental Data Obtained* 



(a) Excess Reactivity Available 

One of the most important parameters required for the high 
power estimation is the maximum excess reactivity available* Since this 
quantity is very small, it cannot be calculated with any degree of ac- 
curacy but may be determined experimentally with reasonable precision from 
the reactor period* The rate at which excess reactivity can be added is 
also of importance* Two forms of reactivity can be accomplished* the 
step input and the ramp input* The ramp input is the usual form of in- 
put* but a step input could be arranged if necessary* 

The method of applying a ramp input would be the insertion of a rod* 
and, since the coarse rod contributes more reactivity per unit length 
than the fine rod* coarse rod movement only is considered* The procedure 
then for determination of excess reactivity and rate of application of ex- 
cess reactivity was as follows: 

With the fine rod fully inserted and the glory hole empty* the re- 
actor period was determined for various positions of the coarse rod from 

1 

the Varian Recorder charts* This experiment was repeated with polye- 
thylene in the glory hole* Finally, the actual rate of coarse rod travel 
was established from the average of a number of runs in which travel for 
a specified stop watch time was measured* The temperature of the water 

1 

Two ion chambers and one BF^ counter are fitted in the reactor for 
monitoring the power level* A Varian Recorder is used to record the 
micro-micro ammeter output for either ion chamber or the count -rate 
meter of the BF^ counter* 
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shield on the day of the experiment was (21 0 2 + 0.5)°C and it is assumed 
that the core temperature was very nearly the same. 

The rate of coarse rod movement was found to be ( 26 o 3 ± 0.2) cm/ 
min. ; this figure coupled with an almost linear relation between insertion 
and excess reactivity provided data for a ramp input. The relation be- 
tween rod insertion and vg- , the ratio of reactivity to delayed neu- 

0 

tron effect iveness & is shown in Figure 3(b). It is noted that the contri- 
bution to reactivity per unit length of the rod is the same for polye- 
thylene in the glory hole as for an empty glory hole for practical pur- 
poses. The linear relationship agrees with the approximately linear por- 
tion of the standard S-curve since the rods of the AGN 201 do not pass 
through the entire core. From the slope of the curve in Figure 3(b) , 
the value of per unit length is approximately 5.61 x 10 ^/cm. 

The function of excess reactivity originates from the equation 



The factor # allows for higher effectiveness of the delayed neutrons as 
compared to the prompt fission neutrons in reaching thermal energy. In 
the experiments for excess reactivity, the periods used were sufficiently 
long so that the first term of the eqxiation was insignificant. Since the 



form as shown in Section 7. Further discussion of the origin and signi- 
ficance of the equation is given in Appendix 1. 

Since the worth per unit length of the coarse rod is approximately 
constant over the last few centimeters of travel, the worth of the 





factor ft* is not well known, 




was computed and can be used in this 
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polyethylene in the glory hole can be related to a length of coarse rod* 
Over a period of time, the coarse rod position with and without polye- 
thylene in the glory hole was noted in connection with the temperature 
coefficient estimation* The difference of positions, the worth of the 
polyethylene is 3.94 ± .03 cm, 

(b) Temperature Coefficient of Reactivity. 

The temperature coefficient experiment was actually a rough 

“4 1 

attempt to check the reported experimental value of (-2,75 x 10 )/°C, 

This experiment assumes that the effect on the core of a constant tempera- 
ture throughout is the same as for the numerically equal average tempera- 
tur^ of a temperature distribution in the core. This is not quite ac- 
curate because a volume element near the centre of the core has a &re&L 
er "importance" than a volume element near the surface of the core and 
is at a higher temperature. If an average temperature is used, it should 
be weighted according to the importance distribution. The experiment then 
serves only to provide an approximation of the magnitude of the tempera- 
ture coefficient of reactivity. 

Over a period of several weeks, the overall temperature of the 
reactor was raised and lowered over as wide a range as possible by ad- 
justing the thermostat of the heating system in the reactor facility. The 
temperature range achieved was only two and a half degrees by this elemen- 
tary procedure, but sufficient difference was noted in the position of 
the coarse rod for criticality to make the plot shown in Figure 3(a) 0 



1 



AGN Report #22, Aerojet General Nucleonics. 
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Temperature (°C) 

Figure 3(a) Temperature Coefficient of Reactivity Determina- 
tion 




Coarse Rod Position Relative to Position for Critical ity (cm) 



Figure 3(b) Excess Reactivity as a Function of Coarse Rod Position. 
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The least square curve through the points has a slope of + ,472 cm/°C 



prisingly well with the value obtained by AGN for a much wider range of 
temperatures o 

At this pointy the questions may be raised as to whether the worth 
per centimeter of the coarse rod determined at 21*2°C is the same for 
other temperatures around 20°C and if the polyethylene has an effect on 
the temperature coefficient estimate as compared to the reactor with the 
glory hole empty* A rerun of the excess reactivity experiment at a dif- 
ferent temperature produced no significantly different results and the 
worth of the polyethylene over the temperature range showed no consist” 
ent tendency to increase or decrease as a function of temperature* 

(c) Drop Test 



the polyethylene in the glory hole was carried out; the primary purpose 
of this experiment was to approximate a negative step input so that the 
computer programs could be checked* This experiment consisted of obtain- 
ing a Varian Recorder chart of the reactor behaviour when the polye- 
thylene was pulled suddenly from the glory hole while the reactor was 
at a steady power* This procedure is similar to the Drop Test method of 
obtaining reactivity for a group of rods (1) (4)« Plotting the logarithm 
of flux versus time over the first few seconds and extrapolating back 
from the approximate straight line slope as shown in Figure 4 provides a 
value of flux 0^ * on the t - 0 axis* Then from the equation 



- H 

corresponding to a reactivity of - 2.75 x 10 /°C. This agrees sur- 



A further experiment providing some insight to the worth of 





1 -H 




17 



t 




18 



Figure 4 Drop Test Simulation! 



a 



I 






| 



a value for of 0*0023 + 0*0002 is obtained which agrees with the 

* 3 " 

value 0*0022 obtained from the product of coarse rod worth per centimeter 
and worth of polyethylene in centimeters of coarse rod* 

(d) Neutron Flux Traverse through the Glory Hole 

The relative thermal neutron flux and the epi-cadmium flux 
through the glory hole are shown in Figure 5* The fluxes were obtained 
by irradiating bare and cadmium- covered indium foils at various equi- 
spaced positions in the glory hole and determining the activity of the 
radiated foils with a scintillation counter* The flux plots were consider- 
ed necessary to obtain an indication of geometric buckling by measuring 
the extrapolated core radius* Also, from the slopes at the interface and 
the curvature in each region, relative values of the diffusion coefficients 
and values of diffusion lengths could be estimated* The experimental pre- 
cision of these flux plots is not particularly good, but the plots serve 
the purposes for which they were intended* 

It is of interest to note that there is a significant dip in the 
thermal flux at the centre of the core due to the thermal fuse so that 
the peak thermal flux is about one centimeter from the center* Also, 
the functions associated with the minor buckling of a two-group analysis 
which are usually positive for thermal flux and negative for fast flux 
at the core-reflector interface are of reversed sign for the AGN 201* 

This fact was useful as a guide to obtaining relative values of dif- 
fusion coefficients and flux equation terms in the parameter determina 
tions discussed in the next section* One further use made of the flux 
plot was a location of the center of the core; first attempts at obtain- 
ing a flux plot indicated that the foil holders in use were off by about 
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Figure 5 Plot of Relative Flux 
in the Glory Hole 
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Distance from Center of Core (\nchos) 



one-half inch in calibration and subsequent corrections were made. 

The data of this section and the data of the next two sections are 
summarized on page 49 0 



21 



5* Theoretical Nuclear Parameters * 



As was mentioned previously in Section 3, two-group theory pro- 
vides a reasonable approximation for reactors moderated by hydrogenous 
materials* If a set of two-group parameters are obtained which predict 
quite accurately several known features of the reactor, it is assumed 
that these parameters and other data calculated from them are good ap- 
proximations for subsequent computations* 

Two-group theory for a reactor with a single reflector region pro- 
vides the following diffusion equations assuming no spatial dependence 
of coefficients within each regions 

D, V 2 <f\(f) - 1‘r^r) 4- v£A(r) = 0 

Core T 

d 2 vV> - i%c f > + =o 

D,V 2 <U?)-E ( ?&(?) — o 

Reflector 

pr^cf) = o 

The solutions of these equations provide equations for the steady-state 
slow and fast flux at all points in the cote and reflector and a critical 
ity determinant from which the critical mass can be calculated* 

(a) Reflector Parameters 

For the reflector, five parameters are required^ the thermal 
and fast diffusion coefficients, the macroscopic thermal cross section, 
the resonance escape probability and the fast neutron removal cross 
section. The removal cross section is obtained from the ratio of the 
fast diffusion coefficient over the ficticious fast diffusion length 
(frequently taken as equal to the Fermi Age)* 



22 



The reflector of the AGN 201 is manufactured from high-density 

3 

reactor-grade graphite with a nominal density of 1 0 7 5 gm/cm e A rough 

3 

check on the density of a sample yielded 1*74 _f *02 gm/cm so that the 

nominal figure is satisfactory for calculat ions* The standard density 

3 

for graphite in most references is 1*60 gm/cm and the tabulated nuclear 
parameters for eommerical graphite are based on this density (1)* Apply- 
ing a correction to the tabulated parameters for density difference* the 
following data are obtained; 



Then. - 

The remaining factor* the resonance escape probability p r * is ef- 
fectively 1.00 since carbon has no resonance peaks in the epithermal 
range (5) and the impurity content of reactor-grade graphite is vety low* 

(b) Gore Parameters - First estimate 

On fortunately $ core parameters are not as readily obtained. 
Polyethylene has been used quite frequently in shielding and as the modera- 
tor of critical assemblies, but* for the most part* the nuclear data for 
water have been used in calculations (with a hydrogen atom density ratio 
correction)* Very little experimental nuclear data on polyethylene it- 
self are available* 

From the official records of "License Material** at the U* S* Naval 

Postgraduate School* the following data pertain to the core: 

235 

Total weight of 0 665*03 gm 



D 



4- 

( 4 ) 



0.838 cm. 



-i 



» 0.000374 cm 

4- a 

D 3 

'J' m 304 cm 
Db ~ ^ ”1 

r 



0.929 cm 
2 



0.00306 cm 
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i 




Percent U 



in uranium 



19.93% 



235 

235 3 

U density in polyethylene 0.055 gm/ cm 

Total weight of components 14781.6 gra 

Total weight of Uranium 3337.67 gms 

The uranium is in the form of uranium dioxide particles compacted with 
the polyethylene. 

235 

The value of 665.0 gms for the weight of U is probably quite ac- 
curate, and for all components, except the thermal fuse, the weight ratio 
235 

of U to total weight is consistent (0.0450 + 0.0003). If the impregna- 

235 

tion density ts 0.055 gm U cc., the overall density is 1.22 gm/cc. 

From the weight of components and the volume calculated from A GN 201 
blueprints, a density of ( 1.22 + 0 . 01 ) gm/cc. is indicated for core 
discs but a slightly lower value is obtained for rod components. Since 
the rods form but a small part of the core and the error involved is less 
than one percent, an overall density of 1.22 gms/cc. will be assumed for 
the polyethylene and uranium dioxide. 

3 3 

The core also contains about 250 cm of aluminum, and about 50 cm 
of pure polyethylene will be in the glory hole to obtain maximum excess 
reactivity. The aluminum has little effect on overall nuclear character- 
istics of the core since both the absorption cross section and scattering 
cross section of aluminum are relatively small. The polyethylene contri- 
butes significantly to the excess reactivity as previously noted, but the 
volume concerned is about 0.47 o of the core volume. The calculations of 
the core characteristics can then be based with little error on the uran- 
ium-impregnated polyethylene with a density of 1.22 gm/cc. On this basis 
the Table 2 data are calculated. Since the volume ratio of UO^ to the poly- 
ethylene moderator is small, the fast and thermal diffusion coefficients 
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Table 2 



Core Material Data 



235 

Weight of U 


665 • 03 gm 


Core material volume 


12,100 cm 3 


Overall density of core 


1.22 gm/cc 


235 

Densities U 


0.055 gm/cc 


238 

U 


0.221 gm/cc 


Uranium 


0.276 gm/cc 


Oxygen 


0.037 gm/cc 


uo 2 


0.313 gm/cc 


Carbon 


0.779 gm/cc 


Hydrogen 


0. 130 gm/cc 


Polyethylene 


0.909 gm/cc 


Volume ratio of UO^ to total 


volume 2.97. 



and the Fermi Age are approximately the same as for Dure polyethy lene c 

Three sources of neutron cross sections are included In the Biblio- 
graphy and are referenced in Table 3 below* Although the 1960 Supple- 
ment to BNL 325 is not included, data from the supplement were used in 

reference (7) * All microscopic cross sections are in barns* 

235 

For U , Westcott 9 s tables (6) will be used since they provide the 
most up to date data on materials with non - — cross sections and aljov 
for the 1/E tail of the neutron di st ribut ion, 

Using the densities of the individual materials and the cross sections 
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as listed* microscopic absorption cross sections averaged over a Max- 
wellian distribution for a temperature of 20°C are obtained from the 



equation 



- 1 f 

l . 1 28 V v/l / 

Table 3 




Neutron Cross sections 





C^cx 






V 


H (gas 


.332(5) (7) 


38(5) (7) 






C 


.0034(5) 


4.2(5) (7) 








.00373(7) 








0 


.0002(5) 


4.2(5) (7) 








.0002(7) 








„ 235 


694(5) 


10(5) 


582(5) 2.07(5) 


2.47(5) 




694(7) 


10(7) 


582(7) 






683.04(6) 




577.01(6) 2.0695(6) 


2.4498(e) 1 


u 238 


2.71 (5) (6) (7) 




.0005(7) 




The cross 


235 

sections for U 


in the above < 


equation are obtained 


f r om We St- 



eott 8 s tables as 



Westcott 8 s values are those to which tabulated cross sections were 
normalized and differ slightly from accepted values*, 
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A, 

o .p ss 562* 5 barns 

A, 

O a « 669*2 barns 



for r s 0,03 and the cat - off. The value of r defines the pro- 
portion of epithermal neutrons in the spectrum; whereas this was not 
measured experimentally for the AGN 201, the value is considered to 
be a reasonable estimate* The S 2 cut - off corresponding to cutting 
off the neutron spectrum above 4*95 KT, is recommended in Westcott°s 
tables for general cases. 

^£ 35 * -1 

Then - 0.0837 cm 



and hydrogen lumped together* 

With these cross sections the product of the thermal utilization 
and the fast neutrons produced by thermal fission per thermal absorp- 
tion in the fuel can be obtained as 



and 



^ 23£T -1 

s 0.0710 cm 

XL q « 0.0231 cm ^ for the oxygen^, carbon 



- 1.594 



From the equation (1) 




1 



Four figure accuracy is not justified but data of the tables main- 
tain the additional figures. 
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I 













and values tabulated with the equation, an approximation can be obtain- 
ed of the fast fission factor for the AGN 201* As expected £ does 
not differ significantly from 1,000. Finally, with the calculation of 



the resonance escape probability " p " using epithermal scattering cross 
sections and Figure 14 of Section 6 - 2 in Reference (1), a value of 



The fast diffusion coefficient for the core material as a function 



where "i" refers to the nuclide. Since . varies with lethargy, a 
numerical integration was carried out using BNL 325 curves and equal 
lethargy intervals. The contribution of oxygen and uranium scattering 
was neglected since the macroscopic scattering cross sections of these 
materials are relatively small. The numerical integration yielded a 
value of 0.623 cm for D^. 

To determine the removal cross section » the effective "fast 

diffusion length" is first obtained. For hydrogenous materials, the cSo- 
called Fermi Age is approximately equal to the square of the "fast 
diffusion length" (2). An experimental value for the Fermi Age of 
polyethylene has not been published but the Age of polyethylene can be 
calculated from that of water because the hydrogen of both accounts for 
most of the macroscopic scattering cross section and the scattering cross 
section of carbon and of oxygen are similar in the epithermal region. 



can be obtained for the core material. 



p c = 0.960 
k ro =*r(fep t = 1-530 



of lethargy may be approximated using the equation 
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The principal difference between polyethylene and water is the atomic 
density* The Fermi Age of a material varies inversely as the product 
of it: density and its averaged microscopic scattering cross section,, 
Roughly then 



r 



2 

s 



and since 



? r, 

2 

'j' = 27.7 cm for water (2) 

= 0.920 for water 
^ = 0.908 for polyethylene 

Z s - l - 50 cm * for water 
27^ ~ 1,78 cm ^ for the core material 
then - 19,8 cm for polyethylene 

The removal cross section is given by 



- d, _ 
r ~ T ~ 



0.0237 



-/ 

cm 



For the thermal diffusion coefficient, it is necessary to estimate 
again from the known data for water. The thermal diffusion coefficient 
for water averaged over a thermal distribution is given as 0,159 cm* 
Correcting for the ratio of atomic densities provides a first estimate 
of 0,136 cm, for polyethylene if the scattering cross section and the 
average cosine of the scattering angles are assumed similar. However, 
the scattering cross section of polyethylene is as much as 15% higher 
than that of water at thermal energies and exhibits more strongly the 
effects of molecular binding. If the discrepency was averaged over a 
thermal distribution covering the entire thermal range, the difference 
might be in the order of 10% resulting in an estimate of 0,123 cm for 
D 2 « As a ch^-ck on this supposition, the relation 
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■ 












3D = -=r- — L 

\ 1 

was averaged numerically over a Maxwellian distribution with kT equal to 
0o0253 ev 4 Linear relations were assumed for (j '^Mo) °f water and poly- 
ethylene and XT t as a function of energy was taken from BNL 325 in the 
case of water and from an unpublished KAPL Paper for polyethy lene. * 

The ratio of the two averaged functions when applied to the diffusion 
coefficient of water (0.159 cm) yielded a corresponding value for poly- 
ethylene of 0ol25 cm 0 This value is assumed as a reasonable estimate 
for purposes of a two-group analysis. 



(c) Checking Parameters in Critical Mass Calculations. 

With the nuclear parameters calculated thus far in this sec- 
tion, a two-group two-region analysis for cylindrical geometry was at- 
tempted for purposes of computing the critical mass as one check on the 
validity of the parameters* The results of the attempt were unsatisfact- 
ory, and subsequent attempts were made with small variations in the para 
meters in an effort to establish whether the poor results were due to 
minor inaccuracies within the tolerances of the basic data. 

The two-group analysis provides a four-by-four determinant, the 
terms of which consist of functions dependent on the geometry and on the 
coefficients of the partial differential equations. This determinant can 
be expanded into an equation w 7 ith the most sensitive function of core 
parameters on one side of the equation and the less sensitive functions 
on the other side (2). 



^Private communication with R. E. Slovacek, Knolls Atomic Power 
Laboratory. 
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Making slight adjustments to various parameters of the core in 
either the determinant or the equation to achieve a balance is very 
laborious and tends to obscure the significance of changes,, The matrix 
formulation of Reference (1) however separates out the parameters of the 
reflector* which are fairly well known* into a reflector matrix* and the 
adjustment of parameters in the core matrix is both easier and more en- 
lightening,, 

When satisfactory results were still not achieved* the analysis 
was repeated for spherical geometry to establish whether the "corners^ 
omitted in the cylindrical geometry were the source of error * These 
corners are shown in Figure 6(a). The cases shown in Figure 6(b) and 
Figure 6(c) can be solved analytically and the case in Figure 6(d) can 
be approached by iteration using the geometries of 6(b) and 6(c) alter- 
nately. However* with a small cylindrical core surrounded by a relative- 
ly thick reflector* the corners omitted in 6(d) are significant. The act- 
ual core dimensions are very close to the optimum case for cylindrical 
geometry* and this in turn has very similar characteristics to a sphere. 

If an equivalent spherical geometry is assumed with the same buckling as 
the cylinder and approximately the same reflector savings* the omission of 
a large part of the reflector in the cylindrical geometry case can be 
avoided. 

However* it was found that no justifiable adjustment of the estimated 
parameters would satisfy the criticality condition. The parameters in 
every case ltd to a core size considerably smaller than the actual core. 

(d) Effect of Spectrum Hardening 

A factor which had not been included in the calculations ( A 
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Figure 6 Geometries for Analytical Solutions 
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is not included in the calculations thus far for comparison purposes) 
was "neutron spectrum hardening"* This additional factor* resolves the 
problem of not predicting tVie correct critical mass* 

Spectrum hardening refers to a shift of the curve for neutron dis- 
tribution as a function of energy towards higher energies,, When the 
moderating material has a constant scattering cross section over the 
thermal energy range, the shift is due to the -~r- or approximately 
absorption cross section of nuclides contained,, In such cases* 
the shift is slight if the absorber concentration is not large,, When 
the scattering cross section changes significantly in the thermal range 

due to effects of molecular binding* the spectrum is further hardened* 
12 3 

Recent papers * 9 indicate that for polyethylene moderated assemblies* 

the shift is very significant and leads to lower averaged macroscopic 

cross sections* Neglecting this effect can result in a large difference 

in the criticality determinant* The spectrum shift in water can be pre- 

4 

dieted quite accurately using Nelkin’s Kernel * but no kernel has been 
developed foi polyethylene* 



D„ R* Bach* S* T Bunch, J* R* Roesser and R* E. Slovacek, Low Energy 
Neutron Spectra Measurements in Polyethylene Moderated Media* Trans- 
actions of American Nuclear Society, Vol* 3, No* 2* Dec* 1960. 

2 

J. R. Beyster* J* L* Wood and p C* Honek, Spatially Dependent Neu- 
tron Spectra, Transactions of American Nuclear Society, Vol* 3, No* 

2, Dec* 1960. 

3 

J* R* Beyster, J* L. Wood, W* M* Lopez and R* B* Walton, Measurements 
of Neutron Spectra in Water, Polyethylene and Zirconium Hydride* 

Nuclear Science and Engineering, Vol* 9* No* 2, Feb* 1961* 

4 

M. Nelkin, Scattering of Slow Neutrons by Water'*, Physics Review* 119, 
1960* 
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A survey of the papers on this subject for polyethylene-moderat- 
ed assemblies indicated that the peak of the distribution should occur 
at 0,029 or 0,030 electron volts for the absorber to scatterer ratio 
of the AGN 201 core material (the unshifted spectrum peak is at 0,0253 
electron volts for 20°C), In calculating macroscopic cross sections 
averaged over this shifted spectrum, an effective neutron temperature of 
336°K to 347°K would be used. These effective neutron temperatures 
correspond to ratios of average neutron velocity to the 2,200 meter/sec 
value of 1,21 to 1,23 as compared to 1,128 for an unshifted spectrum. 

The empirical formula (1) 



0/ 



1. 1 28 + 1.36 



Z»(k.T) 



where qf is the average neutron velocity of the distribution and \T 
is the velocity corresponding to energy 1< 0 T yields a value of 1,23 
which tends to confirm the data taken from spectrum plots. For purposes 
of calculations assume for absorbers 



c! = 



Applying this equation 
of the core materials. 



I- 128 U 



"To j 

T e ff 2200 



To = 293°K 



T 0 



where T is the actual 

I , 22 U T 

temperature 



to recalculation of the macroscopic cross sections 



■^=■735 

.0781 cm 

■?= Z35 

2_ f * .0646 cm 




.0214 cm 



-1 

-1 

-1 



for carbon, oxygen and 
hydrogen combined 
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